Nano- mechanical magnetization reversal 
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The dynamics of the ferromagnetic order parameter in thin magnetic films is strongly affected by 
the magnetomechanical coupling at certain resonance frequencies. By solving the equation of motion 
of the coupled mechanical and magnetic degrees of freedom we show that the magnetic-field induced 
magnetization switching can be strongly accelerated by the lattice and illustrate the possibility of 
magnetization reversal by mechanical actuation. 
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The dynamics of the order parameter of small mag- 
netic clusters and films is a basic problem of condensed 
matter physics with considerable potential for technolog- 
ical apphcations Q. The magnetic field induced rever- 
sal in such systems is perhaps the most active field of 
research. In clusters as small as 1000 atoms the magne- 
tization carries out a coherent motion according to the 
Stoner-Wohlfarth model 01 • In small magnetic wires, 
on the other hand, magnetization reversal is achieved by 
domain walls traversing the sample 01 • With decreasing 
size of magnetic memory cells the fundamental limits to 
the speed and energy dissipation of the magnetization 
switching are important issues. Ingenious mechanisms 
like the so-called precessional switching in which the 
magnetization vector traces straight paths on the unit 
sphere might come close to the optimum in both respects, 
but the magnetic fields cannot be strongly localized, its 
creation therefore wastes energy. An alternative is the 
current-induced spin-transfer torque 0, that can switch 
magnetic layers 0, Q| as well as move domain walls 0|. 
Completely different switching strategies, e.g. using an- 
tiferromagnets 0, attract interest as well. 

Advances in fabrication and detection push the 
fundamental-mode frequencies of nanodevices up to the 
GHz range [l3,EH- Nanomechanics and nanomagnetism 
come together in magnetic resonance force microscopy 
that has already reached single-spin sensitivity It 
has been suggested that mechanical oscillation can be 
used to force coherent motions on uncoupled nuclear 
spins 01 . Here we propose employing the resonant cou- 
pling between magnetic and mechanical degrees of free- 
dom, studied before only in the limit of small magnetiza- 
tion oscillations to accelerate magnetization rever- 
sal and suggest a mechanism to switch magnetization by 
mechanical actuation alone, i.e. without an applied mag- 
netic field. To this end we solve the strongly non-linear 
problem in the limit of weak damping and energy transfer 
analytically. Comparison with numerical simulations in- 
dicates that the solutions are robust beyond their formal 
regime of applicability. 

We consider a small dielectric cantilever with a single- 
domain ferromagnetic layer deposited on its far end (see 




Figure 1: A nano-magneto- mechanical cantilever supporting 
magneto-vibrational modes. On a dielectric substrate (such 
as Si) a single-domain ferromagnetic film is deposited at the 
free end. 

Fig. 1) in the presence of an external field Hq. 

The dynamics of the magnetization M is well described 
by the Landau-Lifshitz-Gilbert equation 0|: 

„ fdM\ 
— = -7MxH., + .Mx(^_j^^^^^, (1) 

where 7 denotes the gyromagnetic ratio. The phe- 
nomenological Gilbert constant is typically a < 0.01, 
and the derivative {^)^^^^ = + ^(_M,x + M.z) is 
taken in the reference system of the cantilever (see Fig. 
1). For small f — fiL), where (p{y) is the torsion an- 
gle at position y of the cantilever, the effective field is 
Hcff — {vMr^ip — vM^)yi -\- vM^Lp'L -I- Hq, where v de- 
scribes the demagnetizing dipolar field [y ~ 47r for our 
geometry) . The coupling originates here from the demag- 
netizing field and since the crystal anisotropy field can be 
small (e.g. in permalloy) it is disregarded initially. 

Without coupling and Gilbert damping the dynam- 
ics of the magnetic subsystem can be solved analytically 
[l^ . leading to the trajectories and oscillation periods 
depicted in Fig. 2. The period of the motion can be 
expressed by elliptic integrals K as 

Ti - 4V2if {v-lv+) I (7V^) ; E < MHo, 

T2 ^ 2V2K (1 - p+ /p_) / (77=pr) ; E>MHo, 
p± = Hi - uE ± HoyjH'l-2vE + v^K'P , 

(2) 

where Emg = —HqM + The strong variation of 

the periodicity has important consequence for the cou- 
pling to the lattice as explained below. 
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Figure 2: The periods Ti (Emg < 1) and T2 (Emg > 1) 
of the magnetization dynamics in Eqs. (EJ in units of 



2tv/ yy\/ {Ho + vM)Hoj as a function of energy in units of 
HqM. The inset shows a plot of typical trajectories at differ- 
ent energies on the unit sphere {vM — IQHo). 



(3) 



The equation of motion of the cantilever is |l 

where C is an elastic constant defined by the shape and 
material of the cantilever (C — ^^da^ for a plate with 
thickness a much smaller than width d and /i is the 
Lame constant), / = j{z'^ + x'^)dzdx ~ ad^/\2 is the 
moment of inertia of the cross-section about its center 
of mass, p the mass density, and (3 a phenomenologi- 
cal damping constant related to the quality factor Q at 
the resonance frequency ti^e as Q = uJe/{2P) (at 1 GHz 
Q ~ 500 01). Note that uje can be also a higher har- 
monic resonance frequency in what follows. The clamp- 
ing boundary condition is p\y=Q = 0. The conservation 
law for the mechanical angular momentum V°'(y) for a 
thin slice at point y £ {0, L} (without magnetic over- 
layer) dV°' {y) /dt — T {y) , where the torque T{y) flow- 
ing into the sHce, is modifled by the coupling to the mag- 
net in a region y G {L, L + AL} (AL is the length of the 
cantilever covered by the magnetic layer) as: 

^"""^^ ' ^ TfieM, (4) 



^V^'(L) + (--)M(L)yj =T(L) 

where Tfleid = M x Hq, V volume of the magnet and 
T\y = —CT{y) is the torque flowing through the can- 
tilever at point y {t = dp{y)/dy). When AL <C L, inter- 
nal strains in the magnetic section may be disregarded. 
The magnetovibrational coupling can then be treated as 
a boundary condition to the mechanical problem [l^ . 
which is expressed as the torque CtI^^l exerted by the 
magnetization on the edge of the cantilever: 



1 ( dm 



(5) 



^ ^ , -7M X Ho 

The effect of the coupling between Eqs. 1^ and © may 
thus be summarized as 



1^1,=, = v^mih^ + 7M X Ho 

Hcff = [vM^p — j/Mx)x + i^M^pz + Ho , 



(6) 



where we introduce the parameter p^ = 
MVy'HQl (vM)l {'y2pILuJe), whose physical mean- 
ing is explained below. In the absence of damping, the 
above system of equations can be obtained as well from 
the free energy: 



F = V{-MiIo+-Ml+i^M^M.ML))- 
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dx. 

We flrst address the coupHng strength of a system 
described by Eq. Q. Consider the two subsystems os- 
cillating at common frequency ui. The total mechanical 
energy is then Eme = pILuj'^p'^, where pq is the maxi- 
mal angle of the torsional motion (that will turn out to 
be identical to the parameter introduced below Eq. Q). 
By equipartition this energy should be of the order of 
the magnetic energy E^g — MVHq. The maximal an- 
gle would correspond to the mechanical motion induced 
by full transfer of the magnetic energy to the lattice. By 
equaHzing those energies we flnd an estimate for the max- 
imal angle of torsion po — ^jMVHo/ {pILuj"^), which at 
resonance is identical to the parameter introduced above. 
The coupling between the subsystems can be measured 
by the distribution of an applied external torque {e.g. 
applied by a magnetic fleld) over the two subsystems. 
The total angular momentum flow into the magnetic 
subsystem by the effective magnetic field is {MV/j)lu, 
whereas that corresponding to the mechanical subsys- 
tem at the same frequency is {pILupo) uj. Their ratio 
is poyvM/H^. The maximum angle po derived above 
is therefore also a measure of the coupHng between the 
magnetic and mechanical subsystems. This estimate is 
consistent with the splitting of polariton modes at reso- 
nance of G = PqU/'^vM/ Hq fI3|. An estimate for a can- 
tilever with p = 2330 kg/m'^ (Si) and d = 1 pm leads 
to pIvM/Hq - M/ {"f'^pd^Ho) 10^3. Decreasing d, p, 
Hq or 1/M is beneficial for the coupling. 

Magnetization reversal by a magnetic field in the cou- 
pling regime can be realized even without any damping 
by transferring magnetic energy into the mechanical sys- 
tem. Since we find that po 'ti 1 for realistic parameters, 
the subsystems undergo many precessions/oscillations 
before the switching is completed. The switching is then 
associated with a slow time scale corresponding to the 
global motion governed by the coupling or a (reintro- 
duced) weak damping relative to a fast time scale charac- 
terized by the Larmor frequency. The equation of motion 
for the slow dynamics (the envelope functions) can be de- 
rived by averaging over the rapid oscillations. To this end 
we substitute Eq. lO , linearized in the small parameters 
a, (3 and (^0, into the equations for the mechanical and 
the magnetic energies: 



(8) 



We focus in the following on the regime Ht^ <C vM , which 



usually holds for thin films and not too strong fields, in 
which the magnetization motion is elliptical with long 
axis in the plane and small even for larger precession 
cones. Disregarding terms containing higher powers of 
Mx and averaging over one period as indicated by (...) 

(^) + (2/3£;„e> = ~Vu (M^M^ip) , 

By adiabatic shaping of time-dependent magnetic fields 
we can keep the two subsystems at resonance at all 
times. The slow dynamics (p(L,t) ~ A(t)e*('^+''/2)* and 
Mx ~ VF(i)e*"* in time domain is then governed by the 
equation: 

ip{L,t) - I(t)e*"* and - W?(i)e*("^-^/2)* corre- 
sponding to 7r/2-shifted harmonics are also solutions, 
and the initial conditions determine the linear combi- 
nation of two envelope functions, i.e. the beating pat- 
tern of two hybridized polariton modes 0|. When ini- 
tially all energy is stored in one degree of freedom, 
is 7r/2 shifted from (piL,t) and A2{t) = 0. Eqs. lHHl 
describe a (damped) harmonic oscillator with frequency 
ujipQ^vMlU^ <^uj when vW"^ <t: MHq. Such os- 
cillatory behavior persists for general angles (except for 
motion with very large angle cones close to the antiparal- 
lel configuration at which the frequency is reduced by up 
to a factor 1/2). This is illustrated by Fig. 3 which shows 
a numerical simulation of Eq. ^ for an undamped sys- 
tem excited at t = by a magnetic field Ho at an angle 
27r/3 with the initial magnetization. The number of pe- 
riods necessary to transfer all energy from one subsystem 
to the other is therefore given by ~ 1/ (^4:ipo ^ vM/Hq^ . 

Eq. ifTTill also shows that for damping constants a> ipq/tt 
OT P/lo > (fio \J vM I Hq I TT the beating is suppressed. 

Fig. 3 illustrates that the mechanical system absorbs 
energy from the magnetic subsystem and gives it back 
repeatedly in terms of violent oscillations that are mod- 
ulated by an envelope function on the time scale derived 
above. When the envelope function vanishes the magne- 
tization is reversed and the systems seems to be at rest. 
However, since the energy is not dissipated, the momen- 
tarily silence is deceptive, and the beating pattern re- 
peats. An efficient coupHng requires that the frequencies 
of the subsystems are close to each other at each con- 
figuration, which was achieved in the simulation by the 
adiabatic modulation of the magnetic field Hq according 
to Fig. 2. However, the reversal process is robust; an es- 
timate from Eqs. (EJ for the necessary proximity of the 
resonant frequencies of the mechanical and the magnetic 

1/2 

subsystems is Aui ~ [vM/Hq) (fio^- In that case the 
above estimates still hold. 
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Figure 3: Time-dependent response of the magnetomechan- 
ical system to an external magnetic field switched on at 
t = 0, in the absence of dissipation. Plotted are the x— and 
z— components of the magnetization [uM = 10-ffo). 




t (units 2ii/ffl) 



Figure 4: Magnetization reversal assisted by the mechanical 
coupling. In the dashed line the external magnetic field is 
reduced at moment t — 3.5 to half of its initial value (at 
a = P = 0). In the dotted and solid lines the field is kept 
constant. The solid line illustrates the effect of a mechanical 
damping of /3 = 0.04 {ijM = lOi^o). 
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Figure 5: Magnetization switching by a pulsed actuation of 
a damped cantilever without external field. Inset shows the 
corresponding magnetization trajectory on the unit sphere. 
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Let us now consider the reversal from the state Mz = 
— 1 by pulsed actuation of the mechanical system with an 
energy less than H^M by twisting the cantilever abruptly 
ait — Q. The mechanical actuation is not essential here, 
but it helps the magnetization to quickly escape the in- 
stable equilibrium point. We can suppress the backflow 
of energy by sufficiently damping the mechanical subsys- 
tem as illustrated in Fig. 4 for /? ~ 0.04 with a vanishing 
Gilbert damping, a — 0. Alternatively we may detune 
the external magnetic field out of the resonance after the 
first reversal. We observe that even without intrinsic 
damping the unwanted "ringing" after the switching is 
strongly suppressed. What happens here is that starting 
with all energy stored in the magnetic degree of freedom, 
the system is brought into resonance only for the time of 
one beating. All energy is then irreversibly transferred 
to the mechanical degree of freedom. 

Finally we propose a non-resonant mechanical rever- 
sal scheme analogous to "precessional" switching Q. The 
effective field Hoff (see Eq. J^J) has a component per- 
pendicular to the plane of the film ^ Mip. Under a 
sudden mechanical twist this component acts like a trans- 
verse magnetic pulse about which the precessing devel- 
ops. The mechanical actuation should be fast, on a scale 
{■jipiyM)~^ , but there are no restrictions on the mechani- 
cal frequency now. We integrate this system numerically 
for a strongly damped cantilever /3/a; ~ 1 and a twist 
of = 0.2. In Fig. 5 we plot the reaction of Mz, rein- 
troducing an easy axis anisotropy described by DM^z, 
{D ~ 0.05, a = 0.01) and in the absence of an external 
field. As in the case of the precessional switching, the 
switching time can be minimized by a more careful ad- 
justment of the mechanical pulse in order to realize the 
optimum "ballistic" path between = ±1. 

Summarizing, we investigated the non-linear dynamics 
of coupled magnetic and mechanical fields for a cantilever 
with a ferromagnetic tip. Conventional magnetic-field in- 
duced reversal can be accelerated by finding new mate- 
rials with higher dissipation, i.e. Gilbert damping con- 
stant a. We propose here three strategies for fast magne- 
tization reversal and suppressed "ringing" using conven- 
tional magnets but opening new dissipation channels by 
i) making use of the additional mechanical damping, ii) 
shaping the external magnetic field pulses, thus quickly 
channeling-off magnetic energy when damping is weak. 
Furthermore, we propose iii) a precessional reversal due 
to a mechanically generated out-of-plane demagnetizing 
field without applied magnetic fields. 

The experimental realization will be a challenge since 
the cantilever has to work at very high frequencies (note 
that coupling to higher harmonics may be a partial so- 
lution) or has to be actuated fast. The magnetic film 



should be small enough to form a single domain. With 
the present pace of progress in nanomechanics and nano- 
magnetism, we are optimistic that these conditions can 
be met in the not too far future. 
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